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ABSTRACT 

Objectives: Depression is a worldwide leading cause 
of morbidity and disability. Genetic studies have 
recently begun to elucidate its molecular aetiology. The 
authors investigated candidate genes of monoamine 
neurotransmission and early environmental risk 
factors for depressiveness in the genetically isolated 
population-based Northern Finland Birth Cohort 1966 
(12 058 live births). 

Design: The authors ascertained and subdivided the 
study sample (n=5225) based on measures of early 
development and of social environment, and examined 
candidate genes of monoamine neurotransmission, 
many of which have shown prior evidence of 
a gene-environment interaction for affective 
disorders, namely SLC6A4, TPH2, COMT, MAOA and 
the dopamine receptor genes DRD1-DRD5. 
Results and conclusion: The authors observed no 
major genetic effects of the analysed variants on 
depressiveness. However, when measures of early 
development and of social environment were 
considered, some evidence of interaction was 
observed. Allelic variants of COMTInteracted with high 
early developmental risk (p=0.005 for rs2239393 and 
p=0.02 for rs4680) so that the association with 
depression was detected only in individuals at high 
developmental risk group (p=0.0046 and (5=0.056 for 
rs5993883-rs2239393-rs4680 risk haplotype CGG 
including Val158), particularly in males (p=0.0053 and 
P=0.083 for the haplotype CGG). Rs4274224 from 
DRD2 interacted with gender (p=0.017) showing 
a significant association with depressiveness in males 
(p=0.0006 and 3=0.0023; p=0.00005 and (3=0.069 
for rs4648318-rs4274224 haplotype GG). The results 
support the role of genes of monoamine 
neurotransmission in the aetiology of depression 
conditional on environmental risk and sex, but not 
direct major effects of monoaminergic genes in this 
unselected population. 
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ARTICLE SUMMARY 



Article tocus 

■ Impact on depression of monoaminergic 
candidate genes with prior evidence of gene- 
environment interaction for affective disorders, 
and of dopamine receptor genes. 

■ Gene-environment and gene-gender interac- 
tions in the aetiology of depression. 

■ Effect of measures of early development and of 
social environment on depression. 

Key messages 

■ Genes of monoamine neurotransmission play 
a role in the aetiology of depression conditional 
on environmental risk, especially in males and in 
individuals at high early developmental risk group; 
in particular, there is evidence of an interaction 
with a C0M7"high-risk haplotype including Val158. 

■ Gender-specific mechanisms and responses to 
environmental effectors are evident in the 
regulation of mood. 



INTRODUCTION 

Depression is a major cause of morbidity 
worldwide, with major depression affecting 
5—7% of the population annually and 16% 
over a lifetime. 1 Although a genetic compo- 
nent in the aetiology of major depression is 
evident with a 40—50% heritability, 2 the 
predisposing genetic background has so far 
remained largely undefined, and recent 
findings from genome-wide association 
studies also point to a complex underlying 
architecture. 3 Depressed patients frequently 
exhibit comorbidities such as anxiety and 
alcohol abuse, 4 and certain personality 
types 5-7 have been associated with depres- 
sion proneness. 
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Strengths and limitations of this study 

■ Depression as defined does not necessarily imply clinical 
diagnosis of major depression but is based on a self-report or 
Hopkins Symptom Check List-25 score. Despite this, the 
prevalence of depressed mood was in the same range as that in 
earlier reports. 

■ There was a notable drop-out rate of about half of the original 
cohort members. 

■ The choice of measures of early development and of social 
environment was limited by the availability of variables 
collected. 

■ Advantages include the availability of longitudinal follow-up 
data starting antenatally, enabling inclusion of the environ- 
mental dimension without recall bias. 

■ The cohort's unique genetic structure with isolation and more 
genetic homogeneity permits identification of genetic-risk loci 
that may be missed when using more heterogeneous 
populations. 

■ The subjects are representative of the population, with all 
cohort members born in the same year and within a geograph- 
ically defined area. 

■ The study sample's size is sufficient for identifying genetic 
variants of moderate impact. 

■ Both genders are represented in almost equal amounts; gender 
differences exist in both depression and temperament traits 
such as harm avoidance. 

■ As the sample is a 1-year birth cohort, genetic effects may be 
isolated from the effects of ageing; some psychiatric traits such 
as harm avoidance are age-dependent. 



Environmental risk factors, in particular stressors 
influencing during development, 8 are considered to 
have a significant impact on the development and course 
of depression. It is likely that many of the genetic risk 
factors for depression interact with the early develop- 
mental environment, but recapture of these interactions 
has remained a challenge for aetiological studies of 
depression. Although the interplay between genes and 
environment has been investigated with respect to 
several psychiatric disorders 9 10 including depression, 
this vast subject still remains largely unexplored. On the 
other hand, addressing the effects of genes and envi- 
ronment on psychiatric morbidity enables us to examine 
the two main constituents in their aetiology. Therefore, 
we wanted to include the environmental dimension in 
our study in order also to explore gene— environment 
interactions (GXE). 

According to the monoamine hypothesis, depression is 
caused by underactivity in brain monoamines, such as 
dopamine, serotonin and norepinephrine. 11 Recent 
results of neuroimaging studies have provided further 
support for this theory. 12 The most solid evidence from 
candidate gene studies has perhaps been obtained for 
the interaction of the SLC6A4 gene for serotonin trans- 
porter and stressful early and current life events, 13 
including positive results from a recent review 14 and 
meta-analysis of all studies to date, 15 although there are 



also contradicting results. 16 Other robust genetic 
findings have been obtained on the COMT gene for 
catechol-O-methyltransferase, an enzyme catabolising 
catecholamines such as dopamine and norepinephrine, 
which has been implicated in depression in conjunction 
with stress, 17 and on the MAOA gene for monoamine 
oxidase A, an enzyme-oxidising neurotransmitter and 
dietary monoamines such as serotonin, norepinephrine 
and dopamine, which has been associated with depres- 
sion in interaction with severity of maltreatment in 
childhood. 17 Furthermore, the TPH2 gene for trypto- 
phan hydroxylase 2, which is the brain-specific form of 
the key enzyme in serotonin synthesis, has been impli- 
cated to interact with stress on disorders of cognitive 
control and emotional regulation, including depres- 
sion. 18 Within the dopamine transmission, the DKD2 
gene for dopamine receptor D2 has been associated with 
depressiveness and anxiety, combined with an effect of 
parenting in childhood, 19 and the DKD4 gene for 
dopamine receptor D4 has been associated with an 
increased risk for obesity in women with seasonal affec- 
tive disorder. 20 Thus, genes from the monoamine 
neurotransmission system are among the most thor- 
oughly studied in psychiatric genetics and in particular 
in the aetiology of mood disorders, and have provided 
perhaps the most robust evidence so far for interaction 
with various types of risk environments, including 
childhood environment. 

We chose to include these candidate genes of mono- 
amine neurotransmission showing prior evidence of 
gene— environment interaction, including SLC6A4, 
TPH2, COMT, MAOA, as well as the dopamine receptor 
genes DRD1—DRD5, in our study on the aetiology of 
depression with a particular focus on their interaction 
with available markers reflecting measures of early 
development and of social environment. The study was 
performed in a sample of 5225 individuals from a large 
Finnish isolated population cohort. As gender is an 
important confounder for depression and at least some 
of the genetic liability is gender-specific, 2 we also exam- 
ined gene— gender interactions in this sample. 

METHODS 
Setting 

We utilised the genetically isolated Northern Finland 
Birth Cohort (NFBC 1966) to investigate the effects of 
candidate genes and environmental risk factors during 
the development on depressiveness. We subdivided the 
study sample based on measures of early development 
arising from the fetal growth environment and neuro- 
logical development during the first year of life 
(measure of early development) as well as from the 
family environment during pregnancy and early child- 
hood (measure of social environment). We examined 
interactions of these measures with candidate genes of 
the monoamine neurotransmitter systems, which have 
prior evidence of gene— environment interactions on 
affective disorders, namely SLC6A4, TPH2, COMT, MAOA 
and the dopamine receptor genes DKD1—DRD5. 
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Study subjects 

The Northern Finland Birth Cohort 1966 (NFBC 1966) 
is a longitudinal 1-year birth cohort from an unselected 
population (N=12 058 live births) comprising inhabi- 
tants of the two northernmost provinces of Finland. 21 
Data collection was begun during the antenatal period, 
and follow-up studies were performed at the ages of 1, 14 
and 31 years. The cohort study was approved by the 
Ethical Committee of Oulu University Faculty of Medi- 
cine, and written informed consent was obtained from 
all participants. 

In 1997, for the 31-year follow-up study 22 all alive 
cohort members with a known address (N=ll 540) were 
sent a postal questionnaire surveying lifestyle, social 
status and health (76% participated), including the 
Hopkins Symptom Check List-25 (HSCL) 23 and items on 
self-reported lifetime depression diagnosis (eg, 'Has 
your doctor ever diagnosed a depressive disorder?'). 
Additionally, cohort members who lived in Northern 
Finland or had moved to the Helsinki area (N=8465) 
were invited to a clinical examination (71% partici- 
pated) with another questionnaire to be filled in later 
and sent to the research group (61% participated). 24 It 
included, among others, a validated Finnish translation 
of Cloninger's Temperament and Character Inventory 
(TCI) questionnaire. 25 

Current depressive symptoms were assessed by the 
HSCL questionnaire, 26 a 25-item shortened version of an 
originally 90-item questionnaire. HSCL contains 13-item 
depression and 10-item anxiety subscales assessing pres- 
ence and intensity of depressive and anxiety symptoms 
during the previous week. Answers are scored on a scale 
from 1 (not bothered) to 4 (extremely bothered). The 
HSCL total score is the sum of items divided by the 
number of items answered. We used mainly HSCL total 
score, as symptoms of depression and anxiety are known 
to overlap significantly. In the post hoc analyses, in order 
to better understand the original association signals, the 
separate HSCL subscales for depressive and anxiety 
symptoms were also taken into consideration. In addi- 
tion to current depressive symptoms (HSCL score) and 
lifetime (diagnosed) depression, we used the TCI 
temperament trait Harm avoidance 5-7 and its subcom- 
ponents as a measure of proneness to depression. 

The subjects (n=5225; 2509 males, 2716 females; 45% 
of the 31 year follow-up study sample or 43% of the 
original study sample) were divided into high- and low- 
risk groups based on the available information reflecting 
measures of early neurodevelopment and of social 
environment (table 1). The markers for the measure of 
high early developmental risk included (1) low birth 
weight (<2500g), 21 considered to reflect suboptimal 
growth environment during fetal life and to increase risk 
for somatic and psychiatric diseases such as depression 
in adulthood; (2) late motor development as reflected 
by first standing later than at the age of 10 months; 28 and 
(3) late development of speech, defined by no words at 
the age of 1 year. 28 If two of these risk indicators were 



present, the subject was classified as having experienced 
a high-risk environment for early brain development. 
The markers for the measure of high social risk envi- 
ronment included the occurrence of two or more of the 
following five indicators for high-risk social environment 
during pregnancy and early childhood: (1) unwanted 
pregnancy (rated by mothers of the cohort members at 
the sixth or seventh month of pregnancy), 29 (2) low 
socio-economic status, shown to be linked with depres- 
sion in the offspring in earlier studies, 30 as defined by 
father's occupation at birth (no occupation, unskilled 
worker, or farmer with area under cultivation under 8 
hectares), (3) single parenthood at birth, (4) low level of 
education of mother (less than 9 years of primary 
school) and (5) low level of information retrieval by the 
mother related to pregnancy and childcare. There was 
no significant drop-out in either of the high-risk groups, 
as 43% and 41% of the individuals of high early devel- 
opmental and social risk groups, and 47% and 46% of 
those of the respective low-risk groups, were available for 
study. 

Genotyping methods 

We investigated genes relevant within the context of the 
monoamine hypothesis of depression: SLC6A4, TPH2, 
COMT, MAOA and the dopamine receptor genes 
DRD1—DRD5 (table 2). The genotyping was performed 
at the Broad Institute (Cambridge, MA) on the 
HumanCNV370-duo chip (Illumina, San Diego, Cali- 
fornia) platform according to the manufacturer's 
instructions. The analysed SNPs included HapMap tag 
SNPs (http://www.hapmap.org/index.html.en) and 
were relatively evenly spaced to cover the genes and 
flanking regions. 

Statistical analysis 

LD structures were determined using HAPLOVIEW. 
Interaction and association/ correlation analyses using 
linear and logistic regression with permutation were 
performed using PLINK Software Package Version 1.04, 
in a stepwise manner to maximise our ability to detect 
associations and to minimise multiple testing. First, 
analyses were primarily performed to identify genetic 
risk variants for current depressive symptoms (HSCL 
score) interacting with measures of early development 
(GXE Dev ) and of social environment (GXEs oc )- For 
variants giving significant evidence of interaction, we 
also performed analyses separately in subgroups at high 
and low risk, respectively. As gender is an important 
confounder for depression, and at least some of the 
genetic liability is gender-specific, 2 we also examined 
gene— gender interactions (GxSex). For variants showing 
significant evidence of gene— gender interaction, we also 
performed analyses separately in males and females. In 
order to achieve a more complete view of the effects of 
the examined genes on depressive symptoms in the 
cohort, we also examined their influence on the gender- 
adjusted HSCL score in the complete sample regardless 
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Table 1 Composition of the study sample from the NFBC 1966 



N 



Hopkins 
Symptom 
Check List 
score>1 .75* 



Measure of early 
development! 



Measure of social 
environment! 



Depression 
diagnosis 



High-risk§ Low-risk nd^j 



High-risk§ Low-risk nd£ 



Males 2509 
Females 2716 
All 5225 



169 (7%) 
269 (10%) 
438 (8%) 



79 (3%) 
136 (5%) 
215 (4%) 



229 (9%) 
193 (7%) 
422 (8%) 



2094 (83%) 
2328 (86%) 
4422 (85%) 



186 (7%) 
195 (7%) 
381 (7%) 



912 (36%) 
1034 (38%) 
1946 (37%) 



1574(63%) 23(0.9%) 
1649(61%) 33(1.2%) 
3223(62%) 56(1.1%) 



*There is prior support for using the Hopkins Symptom Check List score 1.75 as a cut-off when aiming to identify clinical depression. 

t Defined by the presence of two out of three possible indicators for high early developmental risk: low birth weight, late motor development and 

late development of speech. 

t Defined by the presence of two out of five possible indicators for high social risk environment: unwanted pregnancy, low socio-economic 
status, single parenthood, low level of education of mother and low activity for information retrieval by the mother. For further details, see text. 
§Both high early developmental and social risk present in 92 males (3,6%) and 67 females (2,4%). 
INot defined. 



of environmental effectors. Finally, we tested for gene- 
environment correlations (rGE Dev and rGs oc ) and asso- 
ciations of the risk environments with the HSCL score 
(PASW Statistics 18, linear regression model). Second, 
haplotype analyses were performed when two SNPs 
located at physically close vicinity had given association 
signals of p<0.05 when analysed separately. Third, 
genetic variants and haplotypes which had been identi- 
fied in the previous analyses were analysed post hoc with 
respect to HSCL subscales (depressive and anxiety 
symptoms), depression diagnosis and TCI temperament 
Harm avoidance. We report pointwise empirical p values 
generated by PLINK's max(T) permutation (10 000 
permutations) throughout the manuscript, and state 
explicitly where corrected empirical p values are 
reported. SNPs with Hardy— Weinberg Equilibrium 
p values <0.05 were excluded from all analyses. 

RESULTS 

Gene-environment and gene-gender interaction and 
association analyses in relation to the HSCL score 

We examined the effects of nine candidate genes of 
monoamine neurotransmission on current depressive 
symptoms (HSCL score) in a longitudinal population- 
based NFBC 1966 cohort. In particular, we searched for 
evidence of interaction of variants in these genes with 
two measures of early growth, one with indicators for 
potentially disturbed neurobehavioural development 
(measure of early development) and the other with 
risk factors from social environment for normal 
emotional development (measure of social environ- 
ment). The results are presented in table 2 in which 
nominal p values are reported. 

Out of the 69 genetic variants examined, none gave 
a statistically significant association signal with depres- 
siveness or for an interaction with measures of early 
development or of social environment, which would 
survive correction for multiple testing. We observed 
nominal evidence for association with the HSCL score 
(p<0.05) in the complete sample in the cases of 
rsl487275 in TPH2, (p=0.049, 3=0.008), rs4646316 in 
COMT (p=0.026, 3=0.012), rs4274224 and rs4581480 in 



DRD2 (p=0.022, 3=0.011; and p=0.009, 3=0.022, 
respectively), and rsl3106539 in DRD5 (p=0.044, 
3=-0.008). Three variants of COMT and one of DRD3 
showed some evidence of interaction (p<0.05) with high 
early developmental risk with respect to the HSCL score 
(p=0.028 for rs737866, p=0.005 for rs2239393 and 
p=0.020 from rs4680 from COM7; and rs9825563, 
p=0.045 from DRD3). All of these were associated with 
the HSCL score in individuals of the high-risk group 
(p=0.036, 3=0.0414 for rs737866; p=0.008, 3=0.0440 
for rs2239393; p=0.042, 3=0.0320 for rs4680; and 
p=0.022, 3=-0.0396 for rs9825563, respectively). None 
of the variants gave any evidence of interaction with the 
measure of social environment in relation to the HSCL 
score. Five of the genetic variants showed some evidence 
of interaction with gender (p<0.05), including rs737866 
and rs5993883 in COMT and rs4274224 in DRD2. Out of 
these, only rs4274224 was associated at p<0.05 with one 
of the genders (p=0.0006, 3=0.023 in males). The 
evidence for gene— environment correlations (rGE) was 
observed only nominally about rsl906451 from TPH2 
(p=0.035), rs265973 from DRD1 (p=0.047) and 
rs9825563 from DRD3 (p=0.028). Despite a priori 
evidence for the role of the markers which indicate 
a high developmental risk for psychiatric health and well- 
being, namely low birth weight 21 27 and late motor or 
verbal development, 28 there was no correlation between 
these markers and the HSCL score in the present sample 
(p=0.131), whereas the social high-risk environment 
correlated significantly with the score (p=0. 00001). 

Although none of the association findings of these 
primary analyses survived correction for multiple testing, 
post hoc association analyses in gender groups led to 
a finding close to statistical significance, even when 
taking into account the amount of multiple testing 
performed (p=0.0006 for males with rs4274224 in 
DRD2). Furthermore, as there was an accumulation of 
association signals within two highly plausible candidate 
genes, DRD2 and COMT, we proceeded to perform 
haplotype analyses on these genes in order to better 
characterise the allelic variants which yielded the 
observed suggestive associations, and to obtain 



4 



Nyman ES, Sulkava S, Soronen P, et al. BMJ Open 201 1 ;1:e000087. doi:1 0.1 136/bmjopen-201 1-000087 



Early environment, gender and genes in depression 



□ x: 
LU .t; 

CD s= 

i 

. c 

> o 

CD ^ 

Q CO 

X o 

OCO 
^ co 

CD £ 

E 02 

fl 

> T3 
CD C 
T3 CO 

.>- ^-i 



=3 

CO ^ 

S 8 
E 6 

I E 

Co o 

■if 

■| CO 

§£ 

g CO 

il 

c O- 



g *< 

O CD 

C > CO 

- » CT) 

O CO °> 

» £ o 

cj aO 

c cd m 

CD T3 u_ 

O CD CD 



3 ** 

° E 
c o 

O i= 
CD 



O X » 



a co e 
S x « 
S)CD m 

g ? CO 

a a) o 

-O CJ>a) 

CD co E 

w H- o 

5 CD 
O o Si 
■*= CO £ 

« CD r- 

CD ■- 

t 5 CD 

O o ^ 

o CO o 

» 

10 CD .fo 



X 
CD 

c 

o 

o 

CO 

£ 



CM 
Q) 

CO 



CD 

E 
s= 

I 

> 

CD 

'o 
o 

CO 

H — — 

O 2 
"O 
!= 
CO 



CO 

CD 



Q 
LU 

CD 



< 
a. 



X 

a> 
co 
X 
CD 



co 
LU 
X 
CD 



o 

LU 

X 
CD 



o «J 

C CD 



o a 

Q- S3 



a. 

z 

CO 



CD 

£ 
o 

CO 

o 
£ 

o 
— 

o 



CD 

£ 

CO 

c 

CD 

c 

CD 

CD 



CD 

c 

CD 

CD 



I 





^- 




CM 


I s - 




(X) 


co 


o 




CM 


0) 


LO 


CO 


LO 




LO 






CM 


CD 


CD 


co 


CM 


CD 


co 


co 


■* 


CO 


CD 




05 


CM 




a> 


CM 


■sT 


CM 


I s - 




I s - 


CO 


CM 


co 


CO 


00 


CO 


CO 




LO 


CD 


CD 


CO 


LO 


co 


co 


o 


O 


CD 


O 


CO 


CM 


CM 


CM 




I s - 




CM 


co 


co 


co 


CO 


CO 




CD 


LO 


CM 




■st" 


CO 


CM 






co 


© 


O 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 



I s - 


CM 


CD 


to 




LO 


CO 


I s - 


CM 


co 


o 




CO 






CM 


o 


CO 


CD 




CO 


co 


CM 




co 


co 








o 


O 


I s - 


CO 


o 


I s - 


CO 


CD 


co 




CD 




CO 


CD 


I s - 


LO 


CM 


CD 


CM 


co 


CO 




CM 


CO 


CO 




CM 


I s - 




CM 


I s - 


CD 


I s - 




co 


co 


co 


co 


LO 


CD 


co 


CM 


I s - 




CO 


CD 


I s - 


CO 


co 




LO 


LO 


CO 


CO 


CM 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


CO 


o 


co 


o 


CO 


I s - 


CO 


CM 




CD 


I s - 


co 


CO 


* 
* 
« 

O) 


I s - 




CO 


LO 


CM 


co 


o 


o 


co 


-i— 
-i— 

CO 


LO 




co 


o 


CO 


CM 


I s - 


LO 


o 


O) 


CM 


CJ) 




co 


C3) 


o 


I s - 




CO 




o 


LO 


CM 


LO 


CM 


CO 


■st" 


CM 




co 


o 


a) 


CM 


CO 


co 


LO 


CO 


^- 




I s - 


co 


I s -; 




CD 




O 






CO 


I s - 


CO 


CJ) 


(35 


CD 


co 


o 


CM 


CM 


CO 


CO 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 





































000 


























uoo 












uoo 


















OCT) 






uoo 
















■* 


0) 


r- 


LO 


CM 


o 


CM 


CO 


CD 


CM 


CM 


co 


co 


LO 


CD 




CD 


*f 


CM 


o 


LO 


CD 




CO 


CD 






CO 


CO 


CD 


co 


o 


CM 


o 


CO 


CD 




o 


LO 


CO 


co 


CM 


CM 






CM 


CD 


co 


CM 


LO 


o 


CO 


co 


CO 


o 


LO 




I s - 


o 


LO 




^- 




CM 


o 






co 




co 


CO 


CM 


CD 


o 




LO 


o 


CM 


LO 


CD 


o 








O 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 



CO 


CO 


o 


OD 


CM 


r- 


CM 


^1- 


CO 


CD 


T 


co 




■st" 


o 




LO 


CO 


co 


co 


LO 


LO 


CD 


LO 


co 


C35 


LO 


LO 


CD 


CD 


CD 


OD 


CO 


co 


CO 


o 


co 


LO 


CO 


co 




LO 


CD 


CM 


LO 


LO 






C35 


CD 


LO 


co 




o 


CM 


LO 


CO 


CM 






CD 


CD 


CD 




CD 






o 


o 


o 


o 




CM 


CO 


CO 


CM 






CD 




LO 


CO 


CO 


CD 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 



CO 


LO 


CO 


CO 


CO 




I s - 


CO 


CM 


CD 


co 


co 


LO 


CM 




CM 


00 


CO 


co 


LO 


o 


LO 


o 


LO 


T 








o 


CD 


CO 


LO 


I s - 




co 


oo 


CD 


CD 




co 


O 




I s - 


LO 


00 


CM 




CD 


CO 


o 


CM 


o 


oo 


CO 






CD 


CO 




CM 




00 


I s - 


a> 


I s - 


CD 


co 


oo 


CO 


C35 


LO 


CD 


CD 


o 


CO 




CM 


o 


o 


CM 


o 


LO 


LO 




d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 






CD 


co 


LO 


CM 


CM 


CD 


CO 


CO 






oo 




I s - 


LO 


CD 


00 


oo 


C35 


co 


T 


LO 


co 


LO 


o 


o 


CM 






o 








CM 






CM 




CM 




CO 




CO 








CM 


CO 


CO 






CM 


CM 


CM 


o 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


CD 


< 


< 


< 


< 


CD 


< 


CD 


< 


< 


O 


< 


< 


O 


< 


< 


CD 


CD 


CD 


< 


o 


CD 


CD 


< 


< 


CD 


< 


O 


LO 


CJ) 


LO 




CD 


CO 


o 


I s - 




I s - 


i^ 


CD 




CD 








a> 


oo 




co 


oo 




CM 




CO 


LO 


CO 


o 






OD 


CD 






LO 


CO 


CD 


LO 


CM 




LO 


CD 


CO 


CM 


o 


CD 


LO 


CO 


CM 




CO 


CD 




LO 


o 


CD 


a> 


LO 


P-» 


00 


I s - 


CO 


o 


LO 


I s - 


LO 


CO 


I s - 


LO 


I s - 


LO 


O 




CO 


o 


co 




CM 




LO 




CD 


CM 


a> 


LO 


I s - 


LO 


CD 


o 


CO 


I s - 


CD 




co 




co 


co 


oo 


CD 




o 








o 




CM 


CM 


CJ) 




CM 


CO 


LO 


CD 


I s - 


I s - 




CM 


co 








(CD 


CD 


CD 


CD 




o 










CO 


CO 


CO 


CO 


CO 


5- 




LO 


LO 


LO 


LO 


LO 


CO 


CD 


CD 


CD 


CD 


co 


CD 


co 


co 


CD 


h- 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


co 


LO 


LO 


LO 


LO 


LO 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


00 


CO 


oo 


co 


co 


oo 


co 


CO 


00 


co 


co 


oo 


CM 


CM 


CM 


CM 


CM 


I s - 


I s - 


I s - 


I s - 


1^ 




I s - 


I s - 




I s - 































00 




CO 


I s - 


00 


LO 


CD 


CO 


LO 




CD 


CM 


LO 


CO 








LO 




CO 


CO 




co 








LO 


LO 


o 


o 




CD 




I s - 


CJ) 


CD 


CD 


C3) 


I s - 


o 




00 


CM 


O) 


CO 


CM 


LO 


co 


CD 








LO 


O) 


CD 


<*■ 


00 


o 


I s - 


CO 


C0 


LO 


^1" 


CO 


CD 




CM 


CM 


CM 




CO 


LO 


CO 


CO 




co 


CO 




CO 


I s - 




O) 


CO 


CD 




I s - 


CD 






CD 


CD 




O 


co 


1^ 


LO 




CD 


CM 


00 


00 




co 


CO 


CD 


o 


CD 


I s - 


co 




CD 


o 


CD 


o 


CD 


?^ 


CO 


CM 


00 


I s - 


LO 


co 


co 


CM 


oo 


oo 


I s - 




I s - 






CJ) 


CO 


co 


^1- 


LO 


LO 




CO 


O) 


I s - 




o 


o 






CO 




O 


CO 




CO 




CO 


co 


LO 


co 


LO 




CD 


CM 


co 


CD 


CD 


co 


co 


CM 




co 




CM 


00 




CM 




CM 








CD 








CD 


I s - 






LO 


CM 


■st" 








co 


CM 


CO 


CO 


co 


CO 


CO 


co 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


co 


C0 


CO 


co 


CO 


CO 


co 


CO 


CO 


CO 


CO 


CO 



o 



CD 

o 

Q. 

o m 

CD CO 
CO £= 



< 

CD 

o 



CO 



CM 

cr CD 

1 s 

tl 



CM 
X 
0. 



CM 
CM 



CD 
CO 
CO 

CD 

6 1 

o S 

CD £ 

O E 




Nyman ES, 



Sulkava 



S, Soronen P, et al. BMJ Open 201 1 ;1:e000087. doi:1 0.1 1 36/bmjopen-201 1 -000087 



5 



Early environment, gender and genes in depression 
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a maximal amount of information on the nature of the 
associations observed. 

Haplotype analysis of COMT and DRD2 variants in relation 
to the HSCL score 

We performed 2-SNP and 3-SNP haplotype analyses 
combining rs2239393 and rs4680 from COMT and their 
neighbouring variants using the sliding window 
approach. Evidence of association was observed for the 
rs5993883— rs2239393 haplotype CG spanning a region 
from the space between LD blocks 1 and 2 to block 2 of 
COMT (supplementary figure 1) (p=0.0049, (3=0.055), 
for the rs2239393-rs4680 haplotype GG in block 2 
(p=0.0072, (3=0.044) and the rs5993883-rs2239393- 
rs4680 haplotype CGG (p=0.0046, P=0.055) in the high 
early developmental risk group, in agreement with anal- 
yses using single variants (table 3). As rs5993883 from the 
haplotype had also given evidence of interaction with 
gender (table 2), we further examined haplotype associ- 
ation in males and females of the high-risk group sepa- 
rately. We found that the haplotypes increased the risk for 
depressive symptoms in males, but not in females 
(p=0.004, (3=0.083 for rs5993883-rs2239393 haplotype 
CG; p=0.0037, 3=0.072 for rs2239393-rs4680 haplotype 
GG; and p=0.0053, (3=0.083 for rs5993883-rs2239393- 
rs4680 haplotype CGG) (table 3). As is evident from the 
(3-values, each of the haplotypes accounts for more vari- 
ance in depression than any individual constituent SNP. 

Haplotype analysis of rs4274224 and rs4581480 from 
DRD2, which gave evidence suggestive of an association 
with the HSCL score in the complete sample, and of 
their neighbouring variants, gave evidence of an associ- 
ation of rs4648318— rs4274224 haplotype GG spanning 
from block 2 to block 3 of DRD2 (p=0.0007, (3=0.041), 
rs4274224-rs4581480 haplotype GG in block 3 
(p=0.0069, (3=0.022), and rs4581480-rs7131056 haplo- 
type GA spanning from block 3 to block 4 (p=0.0071, 
P=0.022) with the HSCL score. The 3-SNP haplotypes 
rs4648318-rs4274224-rs4581480 haplotype GGG 
(p=0.0027, (3=0.032) and rs4274224-rs4581480- 
rs7l31056 haplotype GGA (p=0.0081, (3=0.021) gave 
evidence of an association in agreement with the find- 
ings from the 2-SNP haplotypes as well as the single 
variants (table 4). As one of the variants contained 
within these haplotypes, namely rs4274224, also gave 
evidence of interaction with gender as well as an associ- 
ation with the HSCL score in males, we also examined 
the association in males alone. The association signal 
became stronger for all of the risk haplotypes, being 
strongest for rs4648318— rs4274224 haplotype GG 
(p=0.00005, (3=0.069). Similarly as for the COMT 
haplotypes, the (3-values imply that each of the DKD2 
haplotypes accounts for more variance in depression 
than any individual constituent SNP. 

Haplotype analysis of COMT and DRD2 variants in relation 
to other neurobehavioural traits 

Encouraged by the findings of the haplotype analyses, we 
tested for associations of haplotypes rs5993883— rs2239393 
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Table 4 Haplotype analysis of DRD2 variants on current depressive symptoms (Hopkins Symptom Check List score) in the 
complete sample from the NFBC 1966 



Males and 
females 



Males 



Gene 


Variant 


Haplotype 


Frequency 


3 


p Value 


P 


p Value 


DRD2 


2-SNP haplotype analysis 
















rs4648318-rs4274224 


GG 


0.05 


0.0409 


0.0007 


0.0694 


0.00005 




rs4274224-rs4581480 


GG 


0.07 


0.0220 


0.0069 


0.0321 


0.0023 






AA 


0.48 


0.0116 


0.0161 


-0.0237 


0.0004 




rs4581480-rs7131056 


GA 


0.07 


0.0220 


0.0071 


0.0322 


0.0026 




3-SNP haplotype analysis 
















rs464831 8-rs4274224-rs4581 480 


GGG 


0.05 


0.0326 


0.0027 


0.0437 


0.0019 




rs4274224-rs4581 480-rs71 31 056 


GGA 


0.07 


0.0215 


0.0081 


0.0317 


0.0033 



Empirical p values based on permutation are reported, with p values <0.05 shown in bold. 



in COMT and rs4648318-rs4274224 in DRD2, as well as 
the single variant rs737866 in COMT with other traits 
related to depression, including the HSCL depression and 
anxiety subscales, depression diagnosis and TCI temper- 
ament trait harm avoidance (table 5). In both genes, it is 
evident that the association with HSCL stems mainly from 
the subscale which reflects symptoms of depression and 
not that reflecting anxiety (with HSCL depression 
subscale, p=0.018, 3=0.075 for COMT haplotype CG and 
p=0.0015, P=0.060 for DRD2 haplotype GG; with HSCL 
anxiety subscale, p=0.288, (3=0.02 and p=0.02 and 
P=0.033, respectively) . We did not detect any evidence of 
an association with depression diagnosis or with Harm 
avoidance or its subcomponents. 

DISCUSSION 

We investigated genetic and environmental risk factors 
for depression in a genetically isolated Finnish birth 
cohort by assessing the relative impacts of mono- 
aminergic candidate genes for depression in groups of 
contrasting (high and low) early developmental and 
social risk. We did not observe any robust genetic effects 
of the analysed variants on depressiveness. However, 
when measures of early development and social envi- 
ronment were considered, some signals for association 
were observed, although none of them survive correc- 
tion for multiple testing. Our study sample provided 
modest evidence of an interaction of COMT with the 
measure of high early developmental risk, particularly in 
males, and a contribution of an allelic variant of DRD2 to 
genetic risk for depressiveness particularly in males 
(table 2). 

The COMT gene encoding for catechol- O-methyl- 
transferase enzyme is among the most investigated genes 
in psychiatric genetics. The enzyme degrades catechol- 
amine neurotransmitters such as dopamine, norepi- 
nephrine and epinephrine by catalysing the transfer of 
a methyl group from S-adenosylmethionine to the cate- 
cholamines. Its enzymatic activity varies according to a G- 
to-A transition at codon 158 in the COMTgene, resulting 
in a valine-to-methionine substitution (Vall58Met) on 
the protein level. 31 The enzyme encoded by the Vall58 



allele has a three- to four-fold higher activity than that 
encoded by the Metl58 allele. Here, we found an asso- 
ciation of the haplotype comprising rs5993883 between 
LD blocks 1 and 2 of COMT, as well as rs2239393 and 
rs4680, which are two variants in virtually complete 
linkage disequilbrium in block 2, with depressive symp- 
toms in high-developmental-risk males (p=0.0053). The 
high-risk haplotype included the high-activity variant 
Vall58 of COMT, the allele G of rs4680. This allele has 
repeatedly been found to be associated with a poor 
response to pharmacological treatment of depres- 
sion, 32 33 and a European multicentre study identified 
an association between that allele and early-onset major 
depression. 34 The Vall58 allele has already been found 
earlier to associate with cognitive deficits including poor 
performance in tasks related to higher-order compo- 
nents of processing 35 and perseverative errors, less effi- 
cient physiological responses in the prefrontal cortex 36 
and even schizophrenia based on a meta-analysis, 37 
although the effect was not significant when studies with 
allele frequencies deviating from the Hardy— Weinberg 
equilibrium were excluded. 

In our study we observed evidence for an interaction 
between COMT and a measure of early developmental 
risk on depressive symptoms. This interaction could not 
be explained through gene— environment correlations. 
Nor were we able to detect a significant correlation of 
the measure of early developmental risk with depressive 
symptoms, despite prior evidence for the role of its 
markers, which were low birth weight 21 27 and late motor 
or verbal development, 28 in decreased psychiatric health 
and well-being, including depression. This finding may 
reflect the presence of other environmental risk indica- 
tors which were not examined in our study. However, 
they may also reflect individual variability in response to 
the risk environment and presence of genetic factors 
(such as the COMT haplotype containing Metl58) that 
may relate to resilience, adaptive changes in regulation 
of emotion reactivity and successful coping with stress. 38 
The observed risk also seemed to arise from an aggre- 
gation of the environmental indicators, as none of the 
risk items separately gave evidence of GXE with the risk 
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variants from COMT or DRD3 (data not shown). This 
could reflect a cumulative nature of these environmental 
influences, such that the effect of one marker may be 
weak, but the accumulated effect of multiple markers, 
together with genetic susceptibility, would be strong 
enough to increase the risk for a deviant development of 
emotional regulation and thus depressiveness. 39 There is 
some prior evidence of interaction of COMT with a risk 
environment on psychosis, antisocial behaviour and 
dissociation. A study on children with ADHD showed 
a gene— environment interaction between the Val/Val 
genotype and low birth weight on early-onset antisocial 
behaviour, 40 and the Vall58 allele was also found to 
associate with cannabis use and psychotic symptoms 41 
and with increasing levels of dissociation in those 
exposed to higher levels of childhood trauma. 42 Inter- 
estingly, a recent report 43 revealed an impact of that 
polymorphism on gender-related patterns of regulation 
of emotions (activation in limbic and paralimbic 
regions) in line with findings of the present study. 

Another main finding of the present study, and statis- 
tically the strongest one, was observed in the dopamine 
receptor D2 gene DKD2, where a haplotype comprising 
the intronic variants rs4648318 in LD block 2 and 
rs4274224 in block 3 was found to associate with depres- 
sive symptoms particularly in males, regardless of their 
early environment (p=0. 00005). Dopamine receptors 
have key roles in a variety of processes in the vertebrate 
central nervous system, and dysfunction in dopaminergic 
neurotransmission may therefore predispose to a variety 
of neuropsychiatric disorders. Among the receptor genes, 
DKD2 has attracted the most attention and has been 
implied to have a role in the aetiology of several psychi- 
atric disorders. However, there are only a few previous 
reports on unipolar depression, including positive, 44 
nominal 45 and negative 46 47 findings, and for results on 
depression conditional on risk environment. 44 46 48 

Our varying results for males and females in general 
imply different mechanisms of mood regulation and 
possible gender-specific responses to environmental 
effectors. Gender differences in depression 2 49 as well as 
in temperament traits 49 have previously been reported in 
various populations, including the current one, 50 and 
the prevalence of depression is higher in women. 51 A 
true gender-specific effect of genetic variants on 
depressiveness would not be surprising, as there is 
evidence of gender differences in dopaminergic func- 
tion 52 that may be oestrogen-dependent. 

It is noteworthy that despite previous reports of the 
5-HTTLPR variant, 13 we did not detect any evidence of 
an association for SLC6A4. Similarly, a recent meta- 
analysis did not find any evidence of an association with 
depression alone, or in interaction with stressful life 
events, 16 although a current review 14 and a meta-analysis 
of all studies to date 15 support the positive association 
findings and the role of 5-HTTLPR and stress in 
depression. The SLC6A4 SNPs included in our study tag 
the 5-HTTLPR well (D >0.9), as determined using 
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genotypes from a population-based Finnish Health 2000 
study. 53 Moreover, the LD measure thus obtained is 
conservative, since in the population under current 
study, LD has been shown to be stronger than in the 
general Finnish population, which was represented by 
the Health 2000 study sample. 54 

We did not use the Bonferroni correction for multiple 
testing, owing to limitations of sample size and expected 
magnitude of gene effects in complex traits. Although 
none of the results from the primary analyses (table 2) 
survive conservative correction, a neurobiological 
a priori hypothesis based on previously published studies 
supports the validity of our most robust findings. It is, 
however, noteworthy that they were observed only when 
the sample was conditioned on measures of early devel- 
opment or of social environment, or gender. Still, the 
strongest association signal, obtained using DKD2 s 
rs4274224 with HSCL score in males (p=0.0006), 
remains close to statistical significance, even when taking 
into account the amount of multiple testing performed. 
The finding was further supported by results of our 
haplotype analysis containing rs4274224, which showed 
a statistically significant association with the HSCL score 
in males (p=0.00005). 

There are some limitations in the present study. First, 
it is notable that depression as defined here did not 
necessarily signify a clinical diagnosis of major depres- 
sion. Instead, it was defined based either on self-report 
or on the score from HSCL, which as a measure has its 
limitations. However, the prevalence of depressed mood 
was in the same range as in earlier reports. 1 55 Second, 
there was a notable drop-out rate among the original 
material of all cohort members. About half of the orig- 
inal cohort members did not participate in this study. 
Finally, when the NFBC 1966 study was initiated, it was 
not possible to predict that an investigation such as the 
present one would one day be conducted. Therefore, we 
are limited by the original choice of variables to be 
collected, and the measures of early development or of 
social environment may only be indicators or markers of 
risk rather than risk factors themselves. 39 It is also 
noteworthy that we did not detect any association of our 
measure of current depression with the measure of high 
early developmental risk, despite it being formulated 
based on previous reports of their effects on psychiatric 
health and well-being. 27 30 However, the effect of 
genetic risk may be modulated by early life stress, even 
though the direct link between early life environment 
and current status would be too weak to be detected in 
our study sample, and this modulating effect may be 
seen in the results of the GXE analysis. 

The current study has several potential advantages, 
such as the availability of longitudinal follow-up data 
starting antenatally, enabling us to include the environ- 
mental dimension without any risk of recall bias. 
Another advantage is the unique genetic structure of our 
study cohort, characterised by isolation, founder effect, 
multiple bottlenecks and more genetic homogeneity 



compared with many other isolates,'' allowing us to 
identify genetic risk loci that may be missed in the 
screening of other more heterogeneous populations. 
Furthermore, the subjects were representative, with all 
cohort members born in the same year and within 
a geographically defined area. 

In addition, the size of the sample is sufficient to 
identify genetic variants of moderate impact. We also 
have both genders represented in almost equal amounts 
(48% males, 52% females), which is notable since 
gender differences are evident both in depression 2 49 
and in temperament traits — for example, harm avoid- 
ance. 49 Furthermore, it is beneficial that the sample is 
a 1-year birth cohort, as it is well established that some 
psychiatric traits, such as harm avoidance 57 of tempera- 
ment, are age-dependent. We can therefore isolate 
genetic effects from the effects of ageing. 

Our results support a modest role of COMT and DKD2, 
two genes of monoamine neurotransmission, in the 
aetiology of depression conditional on environmental 
risk, particularly in males, though not direct effects of 
monoaminergic genes in this unselected population. 
These findings imply that the nature of the role of 
monoaminergic genes in depression should be exam- 
ined further in future studies, and pending replication 
in other, independent population samples. 
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